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Abstract 

The geometrical and electronic structure properties of (100) and (110) silicon nanowires in the 
absence of surface passivation are studied by means of density- functional calculations. As we have 
shown in a recent publication [R. Rurali and N. Lorente, Phys. Rev. Lett. 94, 026805 (2005)] 
the reconstruction of facets can give rise to surface metallic states. In this work, we analyze the 
dependence of geometric and electronic structure features on the size of the wire and on the growth 
direction. 

PACS numbers: 73.22.-f, 81.07.Bc, 81.07.Lk 
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INTRODUCTION 



In recent years one-dimensional quantum systems have attracted great interest as building 
blocks for future molecular electronics devices. Besides proposed applications, there is also a 
fundamental interest in these systems, because as the size is reduced down to the nanoscale 
a description at the quantum level is necessary to account for most of the properties of the 
material. 

Carbon nanotubes are certainly the most studied one- dimensional quantum systems. 
However, the dependence of their electronic properties on their diameter and on their chi- 
rality - and the difficulty of controlling such properties during growth - might make them 
an impractical choice for nanodevice engineering. Semiconductor nanowires seem to obviate 
such an inconvenient and they are becoming the focus of a growing research interest. Sil- 
icon nanowires (SiNWs) are an especially appealing choice, as they offer the best possible 
interface with conventional silicon-based microelectronics. 



SiNWs have been successfully doped with both n- and p-type impurities 



nanoscale logic devices have been demonstrated 
been devoted to the synthesis of SiNW-based 



and 



applications 



sensors 
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as well to device fabrication 



11 



12 



A growing research interest has 
leterostructures - oriented to optoelectronics 
BS- and to the use of SiNWs as™ 
, At present, diameters around 100 nm thick can be almost routi nely 
obtained and a few examples of SiNWs below 10 nm have been published Q Q, Q, Q, 
being 1.5 nm the thinnest wire's diameter reported so far jisl ]. 

Theoretically, at first the attention was directed to the structure of the thinnest possible 



SiNWs. Hollow structures and systems ranging from atomic wires to silicon cluster 
aggregates have been proposed [2^. However, a wealth of experimental results 



15 



soon showed that realistic SiNWs present a well-defined diamond-like bulk core. Since then 



significant efforts have been devoted to the modeling of hydrogen-passivated SiNWs 



2J, 12^, 12211 , i.e. wires where the lateral surface dangling bonds are terminated with 
hydrogens. These are extremely interesting systems for theoretical modeling for different 
reasons: (a) the absence of surface states allows an easier approach to the study of all the 



features related to quantum confinement 



25l l27| . e.g. the energy band-gap is well- 



defined and depends only on the growth direction and on the wire diameter; (b) it is a good 
approach to realistic experimental conditions where different kinds of polluting agents, e.g. 



2 



hydrogen [18[, oxygen 



28| . are present and can react with the surface danghng bonds. 
However, studying the surface physics of these low- dimensional systems without passiva- 
tion is of great interest because of the subtle interplay between an unusually large surface- 
to-bulk ratio and the effect of quantum confinement. In a previous publication j29| we have 
shown that in the absence of passivation, the lateral surface of a (100) SiNW undergoes a 
strong reconstruction which can turn the wire into metallic. The conductive surface states 
form at the {100} facets of the wire, at variance with the infinite Si(lOO) surface which is 
known to be semiconducting. Therefore, the diameter of the SiNW must play a very crucial 
role in the determination of the surface electronic structure. 

In this paper we analyze the diameter dependence of the surface reconstruction and of 
the associated electronic structure (100) SiNWs and we extend our study to wires grown 
along the (110) direction. 

COMPUTATIONAL METHODS 

We have performed density functional theory (DFT) calculations, in the framework of 



the generalized gradient approximation 



orbital 
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GGA) l30||- We have used both numerical atomic 
33l |. combined with norm-conserving pseudopoten- 



32| | and plane-wave methods 
tials of the Troullier-Martins type l2J| and with ultrasoft pseudopotentials, respectively. We 

nnn 

have used a double-C polarized basis set 1^ 1^ and a plane-wave cutoff of 20 Ry for 
the one-electron wave function of the valence electrons. 

All the calculations use periodic boundary conditions, thus the supercell size may restrict 
the number of observed reconstructions. For this reason, to explore a larger set of the 
phase space of the surface reconstruction we have also used a non-orthogonal tight-binding 
model 0, 0] to treat system sizes that cannot efficiently be managed within DPT. Even 
in the cases when the TB surface electronic structure turned out to be only qualitatively 
accurate, the minimum energy geometries proved to be very close to those obtained by DFT. 
In all the calculations discussed, the atomic positions have been relaxed until the maximum 
force was lower than 0.04 eV/A. 

We have studied wires grown along the (100) direction of ~ 8, ~ 11, ~ 15 and ~ 27 A and 
wires grown along the (110) direction of ~ 12, ~ 20, and ~ 26 A, being the ~ 8 and ~ 12 A 
the thinnest possible wires with a bulk core of each growth orientation. For each diameter 
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considered we have carried out calculations with up to four unit cells of the unreconstructed 
wire. The systems with the largest number of atoms (large diameter, many unit cells) have 
been only calculated with TB. The reciprocal space has been sampled with a converged 
number of A;-points, ranging from 1 to 12, depending on the size of the unit cell in direct 
space. 

After converging the atomic positions, we have recalculated the electronic structure with 
a three times thicker mesh of fc-points and an electronic temperature approaching zero, in 
order to obtain a highly accurate band diagram, especially for what concerns the crossings 
of the Fermi level. 



RESULTS AND DISCUSSION 



(100) SiNWs 



The first issue to deal with is the shape of the wire's cross-section. Although we restrict 
to a specif ic g rowth direction, different sections are possible. We have discussed this topic 
elsewhere j40|, concluding that (100) SiNWs favor the formation of {100} facets. On the 
other hand, we found that it is not clear if, at such small diameters, smoothing sharp angles 



increases or not the stability 40|. A discussion of the role of the ed ges can also be found 



\ aiscussion oi tne role oi tne cages c 
I41 I or by Ismail-Beigi and Arias 



in the work by Zhao and Yakobson |4J| or by Ismail-Beigi and Arias j4^. Throughout this 
article, we will focus our attention on wires with smooth angles [see Fig.^(a-c)], even though 
most of the conclusions could be generalized to the case of wires with a square section (see 
Ref . |40| for a complete discussion) . 

The first (100) SiNW that we have analyzed has a diameter of ~ 8 A. This size is well 
below the diameter of the thinnest SiNW grown so far (~ 15 A, see Ref. 13), however, it can 
be a useful case study to elucidate the atomic scale mechanisms that are capable of inducing 
a metallic state at the wire's surface. We have found one only stable reconstruction which 
is illustrated in Fig. |21 It is characterized by one row of Si dimers along the z-axis. The 
formation of Si dimers on the {100} facets is a feature common to (100) SiNWs and is one 
of the characteristic reconstruction patterns of the infinite Si(lOO) surface. 

The band structure diagram, displayed in Fig. |2l presents two states at the Fermi energy. 
It is well-known that DFT, though giving the correct dispersion of bands, underestimates 
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band-gaps. Therefore, especially in the case of marginal crossings of the Fermi level it is 
difficult to guarantee the metallicity of the system. An accurate band-gap estimation can 
be obtained within the GW approximation (see for instance Ref. |27|, where the effect of 
quantum confinement in semiconducting SiNWs is studied). For the purpose of the present 
work, the gap renormalization can be modeled with a scissor operator which rigidly shifts 
upward/downward the states above/below a certain energy. In the present case we cannot 
exclude that the proper scissor operator would open a small gap between the two states that 
in DFT electronic structure crosses the Fermi level, in such a way that the 8 A SiNW would 

hydrogen-passivated wires, we notice that reconstructed wires present a reduction of the 
band gap as compared to the bulk value. 

The most stable reconstruction that we have found for the ~ 11 A wire is shown in 
Fig. |2fa) 13]. It follows the same pattern of the previous wire, featuring one-single row of 
dimers. The differences between these two cases are better understood looking at Fig.HJ^a-c), 
where the unrelaxed geometries are displayed. The thinnest wire can form one only row of 
dimers, involving all of the surface atoms of the {100} facet [panel (a)]. The ~ 11 A SiNW 
can also form one only row of dimers, but it leaves one unpaired row of Si atoms [panel (b)]. 
This has two notable effects: it results in a much more asymmetric in-plane relaxation [see 
Fig. ^b)] and leaves a higher density of dangling bonds. The effects are clearly visible also 
in the side view of Fig. |2l and 01 

The asymmetry in the relaxation originates from the high reactivity of the surface atoms 
of the unpaired row. These atoms cannot form dimers, but their double dangling bonds 
would make the system unstable. The situation is depicted in Fig. HJ the gray atom (in- 
dicated with an arrow) eliminate one dangling bond by breaking a Si-Si subsurface bond 
(red atoms, also indicated with a circle in the picture) and entering it. The golden atoms 
(indicated with a cross in the picture) therein are strictly equivalent to the gray atom, but 
they can combine and form a dimer, according to the conventional mechanism of Si(lOO) 
surfaces. 

The next (100) SiNW that we analyze has already been described in detail elsewhere js^. 
In this case the diameter (~ 15 A) starts to be sufficiently thick to permit to distinguish 
clearly between a bulk region, consisting of atoms tetra-coordinated, and a surface region, 
consisting of atoms with lower coordination: the surface region is the one where the recon- 
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struction occurs, while in the bulk region the atoms maintain the ideal position of the Si 
lattice. 
In Ref . 



29l we have shown that this wire sustains two different competing reconstructions, 
which only differ by 3 meV/atom. While one of them is strongly metallic, with four bands 
crossing the Fermi level, the other is at most semi-metallic. In this case the assignment of 
a metallic character, at least to one of the two stable phases, is robust against failures of 
DFT, because two of the four metallic states are degenerate and thus cannot be shifted by 
a scissor operator without violating correct electron counting 2^. 

This opens up a fundamental question: how can surface states of a {100} facet be metallic, 
while the analogous reconstruction of the Si(lOO) surface is semiconducting? In Fig. El we 
show a top view of the c(4 x 2) reconstruction of the Si(lOO) surface; it can be seen that 
by cutting a strip out of it, one recovers the metallic reconstruction of the {100} facet (the 
one referred to as symmetric in Ref. 12911. However, the coordination of the boundary dimers 
change, because they now miss one of the rest atoms (some of them are indicated with an 
arrow in Fig. E)). This in turn alters the packing density of the surface dimers, especially 
in the case of a nanometric size SiNW where all dimers are at the facet boundary. As 
the diameter of the wire is increased, a growing number of dimer rows will recover the 
coordination that they have on the Si(lOO) surface. Hence, when the fraction of dimers at 
the facet boundary becomes negligible a semiconducting behavior is expected. 

The thicker (100) SiNW that we have studied has a diameter of approximately 27 A. 
Up to this point we have distinguished between wires that have an even or an odd number 
of rows of dangling bonds in the {100} facets, because the different way of reducing the 
number of dangling bonds, i.e. forming dimers or through the mechanisms of Fig. |3J proved 
to be the most prominent factor ruling the reconstruction. In the wires discussed so far the 
dimers formed on the facet were characterized by an altered packing density (the mechanism 
is illustrated schematically in Fig. E}, which ultimately results in the metalization. In this 
wire three row of dimers are formed on each {100} facet, so it is the first among the SiNWs 
that we have considered that features at least one row of dimers which has the same nearest- 
neighbor coordination than the infinite Si(lOO) surface. This situation is revealed by the 
band structure of Fig. IHl where, though still b eing metallic, the band-gap starts opening, as 
can be seen by comparison with the 15 A wire [29]. Another evidence of the size effect can be 
appreciated in the side panels of the same figure. There, we have plotted the wave function 
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of the lower metallic state. It can be seen that the overlap between surface dangling bonds 
decreases for the larger diameter, indicating a higher localization of the surface state. 

Summarizing, the metalization of (100) SiNWs is due to the partial distortion of the Si 
dimers on the {100} facets which alters the overall packing density respect to the infinite 
Si(lOO) surface. When the fraction of dimers with the same coordination of the infinite 
surface increases, i.e when the Si(lOO) packing density is recovered, the SiNW develops a 
semiconducting character. 

(110) SiNWs 

In the case of (110) SiNWs the shape of the wire section is automatically defined by the 
hexagonal symmetry of the [110] plain of bulk Si. The hexagonal shape of the section of (110) 
SiNWs has been also confirmed unambiguously by scanning tunneling microscopy [l^ and 
transmission electron microscopy images. It should be stressed that in this case, contrary 
to (100) SiNWs, the angles between vicinal facets are intrinsically smooth. Hence, the only 
issue concerning the cross-section consists in choosing between two possible configurations, 
shown in Fig. [7| one with rings of hexagons organized around a single hexagonal channel 
and another without the central hexagon. The tests that we have performed indicated that 
the arrangement in Fig. U][h) is favored. Therefore, this is the pattern that we have adopted 
to generate SiNWs with growing diameters (Fig. |HI). In these SiNWs the unrelaxed lateral 
surface is made up of {111} and {110} facets. Therefore, the competition mechanism that 
rule the reconstruction is expected to be markedly different to the case of (100) wires. 

The thinnest (110) SiNW that we have studied turns out to have two competing re- 
constructions, which are shown in Fig. El After the relaxation it is hardly possible to 
still distinguish facets with a well-defined crystallographic orientation [see the top view of 
Fig. Efd)]. The difference between to the reconstructions found emerges at the {lll}-like 
facets, where the periodicity of the protruding atoms change, as it can be observed in the 
space-fill representation of Fig. El The difference in energy is negligible, so that both recon- 
structions are expected to occur with the same probability. Most important, both phases 
are semiconducting with a small band-gap, approximately one third of its bulk value. 
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The ~ 20 A has similar features. It presents two reconstructions with a level of stability 
of the same order (the difference is around 3 meV/atom). In this case too, the rows of pro- 
truding atoms can arrange in different ways: one single row in the case of the reconstruction 
with the smallest motif; a second row which alternates protruding to non-protruding atoms 
in the other reconstruction. Contrary to some of the cases described in Sec. „ the differ- 
ences in the reconstruction alter only marginally the band structure and a semiconducting 
character can be safely attributed to both phases. 

Curiously, it turns out that the number of protruding atoms is always a constant. In the 
smallest reconstruction [Fig. fToT a)] none of the atoms protrude outward, while all of them 
do in the opposite facet (not shown here). The opposite point of view was chosen for the 
unit-cell reconstruction of Fig.lHfa), where all the atoms are protruding outward, while none 
of them does in the opposite facet. In the same way, the more varied motif of the extended 
reconstructions of Fig.Efb) and Fig. [TUfb) is echoed in the opposite facet, where a depressed 
atom corresponds to a protruding atom and viceversa. 

The largest (110) SiNW that we have analyzed has a diameter of ~ 26 A. A significant 
difference with the thinner wires described here, is that the diameter has reached a critical 
size so that not only {111}, but also {110} facets are present. The latter reconstructs 
forming a long trough all along the wire axis, while the {111} have a patter similar to those 
seen so far, with a row of protruding atoms. 

A common feature that clearly emerges is that, although there is a dependence of the 
band-gap size on the wire diameter, in none of the (110) SiNWs that we have studied 
the facet reconstruction gives rise to surface metallic states like in the case of wires grown 
along the (100) direction. From the discussion of the results it appears that these wires are 
dominated by {111} facets and {110} facets form only on thick wires. 



The found geometries - and their semiconducting character - a gree with t 



le TT-bonded 



chain model typical of the 2x1 reconstruction of the Si(lll) surface |45|, |46|, |47|, |48| . The un- 
reconstructed surface (the facet in the case of the wire) features one dangling bond per atom 
and the system spontaneously reduces its symmetry to lower the energy of the half-occupied 
states. This is possible because the surface dangling bonds reside on nearest-neighbor atoms 
and the bond-orbitals are close enough to interact significantly stabilizing the surface and 
introducing a band-gap between occupied and unoccupied states Other reconstructions, 
notably the 7x7, are known to be metallic due to the presence of adatoms that have not 
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been considered in the present study. Adatoms could lead to metallic - and possibly more 
stable - phases also in the case of nanowire facets and their influence on the electronic and 
geometrical structure of the reconstruction should be addressed. However, under a certain 
diameter the facets are not large enough to host such extended reconstructions, but the 
existence of adatom phases cannot be ruled out. 

Further insight on the mechanisms that rule the surface physics of these wires can be 
obtained from the analysis of Fig. ^2 In the left panel we have plotted the energy gained 
with the relaxation: 

e,ei = {Elll - EZn/n, (1) 

where E'^^Jr and E^^^^'^ are the total energies of the reconstructed and unreconstructed geom- 
etry (shown in Figs. Q and |H1) of a wire of n atoms. The first thing to observe is that with 
the relaxation (110) wires lower their energy less than (100). This is in agreement with the 
reconstruction mechanisms previously discussed: (110) wires simply rearrange the topology 
of the surface dangling bonds to allow the formation of a 7r-bonded chain; (100) wires, on the 
other hand, go through a much stronger reconstruction of the facet, involving the formation 
of new bonds. 

Both curves are fairly linear, with the remarkable exception of the 11 A (100) SiNW. 
In the above discussion (see Sec. ^ we have suggested that facets of (100) SiNWs with an 
even number of dangling bond rows should be favored. In the opposite case (odd number 
of dangling bond rows, like the 11 A wire), the system has to spend a larger amount of 
structural energy to readjust the facet geometry, because it cannot form ordered rows of 
dimers. This intuitive idea is confirmed by Fig. ITTT a) and becomes clearer by defining the 
surface energy as: 

Csurf = (-E'tot — '^bulkebulk)/'T-surf, (2) 

where ebuik is the energy per atom of bulk Si and nbuik and ngurf the number of bulk and 
surface atoms, respectively [4^. In this way, in the graph of Fig. ITlTb) we isolated pure 
surface contributions and the anomaly of the 11 A wire emerges in a more evident way. 
As a concluding remark, we discuss the cohesion energies defined as 43^ : 

ecoh = (^tot - nesi)/n, (3) 

where Etot is the total energy of the wire, n is the number of atoms and esi the energy of 
the isolated Si atom. The cohesion energy j43| , being directly derived from the total energy, 
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allows to compare the relative stabilities of the systems studied and which structures will 
be favored, in absence of other constraints, in the growth process. The smoother angle 
between vicinal facets favors (110) growth orientation to the sharper angle of (100) SiNWs. 
When the diameter increases the difference tends to disappear: (100) wires reach a critical 
cross-section that allows to form smoother transition facets, while thick (110) wires start to 
develop their own {100} facets, which are expected to contribute with a surface energy very 
similar to the {100} facets of the (110) SiNWs. 

Summarizing, the (110) SiNWs that we have analyzed turned out to be always semi- 
conducting. The danghng bonds on the {111} facets cannot give rise to the formation of 
dimers because, due to the different symmetry of the cleavage, they are second neighbors. 
Incidentally, if they could form dimers this would result in a complete passivation of the wire 
because each atom on the facet surface has only one dangling bond per atom. However, the 
dangling bonds are close enough to interact and to form a 7r-bonded chain typical of Si(lll) 
2x1 surfaces. 

CONCLUSIONS 

We have performed first-principles electronic structure calculations to study the surface 
reconstructions of (100) and (110) of different diameters. 

SiNWs grown along the (110) axis are dominated by {100} facets that exhibit the typical 
reconstruction pattern of the Si(lOO) surface, consisting in the formation of rows of buckled 
dimers. In the case of the thinner wires, the altered packing density of the facet dimers 
gives rise to delocalized surface states that result in a metalization of the surface. Increasing 
the diameter, the coordination of the infinite Si(lOO) is recovered. The largest (100) SiNW 
studied, though still metallic, gives a clear indication of a band-gap opening. 

The case of SiNWs grown along the (110) axis is different. Thanks to a different cross- 
section shape, the neighboring facets can match through a smooth angle and they arc not 
significantly distorted with respect to their homologous infinite surfaces. Hence, the {111} 
facets that dominate the reconstruction follow the 7r-bonded chain model of Si(lll) surfaces, 
conferring to the wire a semiconducting character. The band-gaps are found to be smaller 
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than the bulk value as an effect of the formation of the 7r-surface states, while as it is well- 
known, H-passivated wire of similar sizes result in widened gaps as a consequence of quantum 
confinement. Analysis of the energetics of the wires indicates that at small diameters the 
growth of (110) SiNWs is favored over (100) SiNWs. 
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y Ciencia. Computational resources at the Centre Informatique National de I'Enseignement 
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Orientation diameter [A] Ecoh [eV] 



(100) 



11 
15 

27 



-3.75 
-3.86 
-3.99 

-4.16 



(110) 



12 
20 
26 



-3.93 
-4.08 
-4.14 



TABLE I: Cohesion energies of the SiNWs studied. For thin diameters (110) are favored, but as 
the size is increased the difference tends to disappear. 
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FIG. 1: (Color online) Section views of the unrelaxed [(a) to (c)] and relaxed [(d) to (f)] (100) 
SiNWs with a diameter of ~ 8 ~ 11 and ~ 15 A. The facet geometry is determined by the number 
of rows of dimers (indicated by the brackets) that can form: (a,d) one-single row of dimers can form; 
(b,e) the formation of a second row of dimers is frustrated and the presence of an unpaired row of 

atoms results in an asymmetric relaxation of the facet; (c,f) two rows of dimers form, yielding a 
trough in the middle of the {100} facet. 
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FIG. 2: (Color online) (100) SiNWs of ~ 8 A diameter: side view of the reconstructed (pseudo) 
{100} facet. 
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FIG. 3: (Color online) (100) SiNWs of - 11 A diameter: side view of the reconstructed (pseudo) 
{100} facet. 
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FIG. 4: (Color online) The mechanism that rules the facet relaxation when not all the dangling 
bonds can form dimers. The gray atom (indicated with an arrow) is unpaired and cannot form a 
dimer; to reduce the number of dangling bonds breaks a bond of the underlying layer {red atoms 
(indicated with a circle) and enters it. The golden atoms (indicated with a cross) form a surface 
dimer according to the conventional mechanism of Si(100)-like surfaces. 
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FIG. 5: Top view of the Si(lOO) c(4 x 2). White and gray balls represent the (asymmetric) top 
dimers, dark balls are the second layer atoms. A strip equivalent to the {100} facets of the ~ 15 A 
(100) SiNW is drawn. Each of the atoms represented by white balls have only one second layer 
atom where to rest onto (the missing ones are indicated with an arrow). This induces a tilt of 
the dimer axis when it relaxes, thus achieving a different packing density than the infinite Si(lOO) 
surface. 
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29| and (b) the 



FIG. 6: (Color online) Band structure diagram of (a) the 15 A SiNW of Ref. 
27 A SiNW. The comparison between the two panels shows the trend of band-gap opening as the 
diameter is increased. 
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FIG. 7: (Color online) Possible geometries considered for (110) SiNWs. The section that features 
a central hexagon [shown in panel(b)] turned out to be the most stable. 
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FIG. 8: (Color online) Section views of the unrelaxed [(a) to (c)] and relaxed [(d) to (f)] (110) 
SiNWs with a diameter of ~ 12 ~ 20 and ~ 26 A 
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(a) (b) 

FIG. 9: (Color online) {111} facets of the two stable reconstructions that we have found for the 
12 A (110) SiNW. The different unit cells of the facet are indicated. The reconstructions have 
extension of (a) c and (b) 3c, where c is the axial lattice parameter of the SiNW. In (a) the 
reconstruction features one row of aligned atoms, while in (b) only one atom each three (indicated 
with an arrow) protrudes outward. 
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(a) (b) 

FIG. 10: (Color online) {111} facets of the two stable reconstructions of the 20 A (110) SiNW. 
The different unit cells of the facet are indicated. The reconstructions have extension of (a) c and 
(b) 2c, where c is the axial lattice parameter of the SiNW. In (a) the reconstruction features one 
row of aligned atoms, while in (b) only one atom each two (indicated with an arrow) protrudes 
outward. 
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diameter [Ang] diameter [Ang] 



FIG. 11: (Color online) Energy gained with the relaxation, Eq. ^ (left panel) and surface energy, 
as defined in Eq. |2l (right panel) . 
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